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. A simulation study of the energy released in air due to the development of an 

extensive air shower has been carried out usine the CORSIKA code. The contribu- 
tions to the energy release from different particle species and energies as well as the 
■ typical particle densities are investigated. Special care is taken of particles falling 

^ I below the energy threshold of the simulation which contribute about 10% to the 

total energy deposition. The dominant contribution to the total deposition stems 
from electrons and positrons from sub-MeV up to a few hundred MeV, with typical 
transverse distances between particles exceeding 1 mm for 10 EeV showers. 
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P ■ 1 Introduction 

During the shower process initiated by a primary cosmic ray in the atmosphere, 
^ . in general, only a small fraction of the initial energy reaches the ground as 

high-energy secondary particles. Instead, most of the primary energy is re- 
leased in the atmosphere by ionization and excitation of the air molecules. A 
tiny fraction of order 10~^ is emitted as fluorescence light. As shown by the 
Fly's Eye experiment [1] and its successor, the High Resolution Fly's Eye [2], 
the fluorescence component of an air shower can be used to detect cosmic 
rays at energies exceeding about 10^^ eV. Also the FD detectors of the Pierre 
Auger Observatory [3] and planned experiments such as EUSO [4], OWL [5], 
and Telescope Array [6] are based on this detection technique. In contrast 
to measurements with particle detectors on ground, the observation of the 
longitudinal shower curve in the fluorescence light provides calorimetric infor- 
mation and therefore conclusions about the primary energy which are largely 
independent of the primary particle type and of unknown details of hadronic 
interactions at these extreme energies. 
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In the classical approach of reconstructing the primary energy from the mea- 
sured data [7,8], the amoimt of light emitted by the shower at a particular at- 
mospheric depth X (determined from the observables by "backtracing" based 
on detector calibration, atmospheric corrections and previously performed ge- 
ometrical reconstruction) is converted to a number of charged particles Nch 
assuming a mean ionization rate. A dependence of the fluorescence yield, i.e. 
the fraction of the released energy that is emitted as fluorescence light, on tem- 
perature and density is taken into account according to [9]. Finally, Nch{X) 
is integrated over the path-length and multiplied by a mean ionization loss 
rate and a correction for missing energy is made. Air shower simulations are 
invoked to obtain average values for the mean ionization rate and for the 
correction of the missing energy [7]. 

For such primary energy determinations, it is assumed that the amount of 
fluorescence photons produced locally per unit length dNfi / dl is proportional 
to the local energy release dEj-ei/dX in air. 



with Pair{h) — dX/dl being the atmospheric density, where X and I are mea- 
sured along the shower axis. The fluorescence yield y{T,p) of air (in units 
of emitted photons per released energy) is a key quantity for evaluating the 
calorimetric shower energy. It has to be determined by laboratory measure- 
ments with individual particles as projectiles [9,10]. Due to quenching effects, 
y is dependent on environmental conditions. Temperature and density have 
to be varied according to realistic atmospheric conditions with air as target 
material, maybe even taking humidity effects into account. A considerable ex- 
perimental effort is started to determine y(T,p) more precisely [11]. Different 
experiments, some of them at accelerator facilities, with various measuring 
conditions are being planned. Especially, the assumption of a proportionality 
between fluorescence light production and deposited ionization energy, which 
has been justified to some extent by previous fluorescence yield measurements 
[9,10] will be checked thoroughly. 

The role of the total deposited ionization energy for the calculation of flu- 
orescence light makes a precise knowledge of the dominant energy deposit 
processes in air showers important. In Monte Carlo simulations one has to 
apply a low-energy cutoff for explicit particle tracking. Here we introduce a 
method for calculating the total ionization energy loss independently of the 
particular, applied simulation cutoff. 

In this work we investigate the characteristics of the energy release in air show- 
ers. For this purpose, the CORSIKA simulation code [12] has been adapted to 
quantify the contributions of different particle species in an air shower to the 
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energy release. The electromagnetic component as the dominating one will be 
examined in more detail. 

Our results are also of importance for the planning and interpretation of labo- 
ratory fluorescence yield experiments, as we determine the energy ranges that 
dominate the ionization energy deposit in air showers. 

The plan of the paper is as follows. In Chapter 2 it is described how the energy 
release is calculated in CORSIKA. This includes the definition of a releasable 
energy for particles which are discarded in the Monte Carlo process since 
their kinetic energy is below the simulation energy threshold. In Chapter 3, 
the longitudinal profile of the total energy release and the contributions of 
different particle species are discussed. The lateral energy deposit distribution 
together with the densities of the shower particles are analyzed in Chapter 4. 
The energy spectra of electrons and positrons and their importance for the 
total energy release are investigated in Chapter 5. 



2 Calculation of the energy release with CORSIKA 

2.1 The air shower simulation program CORSIKA 

The Monte Carlo program package CORSIKA [12] is designed to simulate 
the development of extensive air showers induced by various types of primary 
particles (photons, hadrons, nuclei ... ) in a wide energy range up to the 
highest energies. For the particle interactions, external state-of-the-art codes 
are employed. Electromagnetic interactions are simulated using an adapted 
version [13] of the EGS4 code [14] , which includes the Landau- Pomeranchuk- 
Migdal effect [15]. 

To describe hadronic interactions at the energies relevant in this paper, ac- 
celerator data have to be extrapolated by several orders of magnitude in en- 
ergy and into a forward kinematic range unobserved in collider experiments. 
Therefore, CORSIKA offers a choice of various hadronic interaction models 
which differ in predictions of certain air shower characteristics. For the present 
investigation, however, the uncertainty related to the modelling of hadronic 
interactions turns out to be of minor infiuence, as the main characteristics of 
the energy release in showers are determined by electromagnetic interactions. 
For the calculations presented in the following, the model QGSJET 01 [16,17] 
has been employed for hadronic interactions with energies Eiab > 80 GeV, 
while the GHEISHA routines [18] have been used to treat hadronic collisions 
at lower energies. 
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To reduce the computational effort in CPU time, the technique of particle 
thinning [19], including weight limitation [20,21], is applied. More specifically, 
a thinning level of 10~^ has been chosen, i.e. no thinning occurs for parti- 
cles with energies exceeding lO^^E'o, with Eq being the primary energy. In 
case of thinning, the bulk of secondary particles produced in an interaction 
is discarded, and only "representative" particles are sampled. A weight factor 
is assigned to the selected particles to keep energy conservation [19]. Upper 
weight limits of 10~^£^o/GeV for electromagnetic particles and 10~^£^o/GeV 
for muons and hadrons have been chosen to keep artificial fluctuations intro- 
duced by individual particles sufficiently small [22] for the analyses presented 
in this paper. 

In CORSIKA simulations, an (adjustable) energy threshold Ethr is adopted 
for the shower calculation, i.e. particles are followed explicitely for E > Et^^ 
and discarded for smaller energies. This is due to the fact that a detailed cal- 
culation down to smallest energies, for instance to the eV-range, seems both 
hardly possible (CPU time) and hardly necessary for "classical" air shower 
experiments measuring surviving particles on ground. In the realm of fluores- 
cence light, however, low-energy particles also contribute to the ionization and 
excitation of air molecules. Therefore, two categories of shower particles are 
distinguished for calculating the energy release: 

• particles above simulation threshold that are tracked in detail (section 2.2), 

• particles below the simulation threshold that are discarded from explicit 
tracking in the further simulation process (section 2.3). 

The final energy release dE^ei/dX is the sum of these two contributions, which 
will be discussed in the following. 

2.2 Ionization by explicitely tracked charged particles 

The continuous ionization energy loss dEi/dx of a single charged hadron or 
muon traversing matter of thickness dx along its track is calculated by the 
Bethe-Bloch stopping power formula 



where (5 = v/c\s the velocity of the particle in the laboratory in units of the 
velocity of light, 7 is its Lorcntz factor, and z is the charge of the ionizing 
particle in units of e. The two constants Ki = 0.153287 MeV g^-'^cm^ and K2 = 
9.386417 are derived from the tables [23] for dry air. Additionally the EGS4 
routines modified for CORSIKA take into account the pressure dependent 
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Sternheimer correction [24]. 



The transport of electrons and positrons witli ionization energy loss is treated 
within CORSIKA in great detail by the EGS4 routines [14] . The contributions 
of the continuous energy loss stem from the soft bremsstrahlung of photons 
below the simulation threshold and from sub-threshold energy transfer, in- 
cluding ionization energy loss, to atomic electrons. For the latter contribution, 
EGS4 uses the formulae recommended by Berger and Seltzer [25] to apply the 
concept of "restricted stopping power" [26]. 

2.3 Treatment of particles below the simulation threshold 

The amount of fluorescence light produced by sub-threshold particles depends 
on the reactions the specific particle would suffer on its further way, and there- 
fore on the particle type. Antiparticles, for instance, will annihilate and may 
finally release a much larger amount of energy than the (usually quite low) 
kinetic energy. Therefore, it is useful to define, similar to the concept of re- 
stricted stopping power, the quantity releasable energy Er for each particle 
species. This releasable energy consists at least of the kinetic energy of the 
particle. Its contribution to the total energy release is added during the sim- 
ulation process at the position in air where the particle is discarded. 

Stable particles such as electrons cannot release more than their kinetic energy 
{Er = Ekin)- Due to annihilation, positrons have a larger releasable energy 
{Er — Ekin + 2-me — Ekin + 1-022 MeV). However, if the resulting pair of an- 
nihilation photons is above the photon energy threshold, no energy deposition 
takes place but the photons are treated explicitely in the further simulation 
process. 

Muons and mesons are unstable and can release also part of their rest mass into 
ionization. Some part will be carried away by the decay neutrinos, however. 
Antibaryons arc assumed to annihilate with a nucleon under emission of several 
pions, so their releasable energy, apart from Ekin, is increased by the anti- 
baryon rest mass and nucleon rest mass. In case of unstable particles and the 
antibaryons, we assumed an effective fraction of about 1/3 of the releasable 
energy to be taken into account for the total energy release, and the remaining 
2/3 to be "lost" mainly in the neutrino channel. The simplified treatment of 
these particle species below the simulation energy threshold seems justified as 
only the contribution of electrons and positrons below threshold turns out to 
be significant for the present analysis. 

The assumption of locality of the energy deposited by the sub-threshold par- 
ticles is valid for conventional threshold settings in the simulation [22]. While, 
for instance, a usual choice of the electron energy threshold is about 0.1 MeV, 
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Fig. 1. Longitudinal shower development: Energy release of electrons plus positrons 
and of muons plus hadrons. Additionally, the contribution of electromagnetic par- 
ticles below the simulation energy threshold Efhr = 0.1 MeV is given. 

even for relatively large threshold values in the MeV range the particles would 
be stopped within a depth of only a few g cm~^. 



3 Longitudinal development 



Air showers initiated by different primary particle types (proton, iron nu- 
clei, and photons) have been calculated with various primary energies (10^^ — 
10^° eV) and shower zenith angles. The main conclusions discussed in the fol- 
lowing, however, are largely independent of the primary parameter choice and 
also of the shower-to-shower fluctuations (see also [27,28]). The distributions 
given in the Figures illustrate the case of an iron-induced event of primary 
energy 10^^ eV with 45° inclination. 

The longitudinal development of the energy release dErei/dX of the shower 
in air is presented in Figure 1. The definition of the path length dX deserves 
explanation. Plotted is the energy release in the layer of air between slant 
depths X and X + dX, i.e. in the direction of the shower axis. A particle 
propagating with a non-zero angle towards the shower axis thus travels 
through an effective amount of matter of dx = dX/ cos('(9) while traversing 
this slant depth interval. The corresponding larger energy loss, or equivalently 
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the increased total tracklength of the shower particles due to their angular 
spread, is taken into account in CORSIKA. This effect has also been pointed 
out recently in [29]. 

In Figure 1, different contributions to the energy release are displayed. As 
expected, the main contribution stems from electrons and positrons, the most 
numerous charged particles. Around shower maximum, less than 2-3% are 
resulting from muons and hadrons. Thus, electromagnetic particles should be 
the main target for the study of energy release. 

A fraction of about 10% to the energy release is due to discarded electrons 
and positrons below a simulation energy threshold of 0.1 MeV. This value 
is in agreement with findings in [7]. High-energy electrons produce via M0ller 
scattering a significant number of low-energy electrons. Additionally, the many 
low-energy shower photons transfer a considerable energy fraction by Compton 
scattering to low-energy electrons. 

A closer inspection shows that the development of the hadronic and muonic 
parts differ somewhat from the electromagnetic ones: At very early (and also 
late) development stages the electromagnetic fraction to the local energy re- 
lease is decreased. For the current analysis, this is of minor interest, however, 
as the fluorescence light production at these shower stages is marginal com- 
pared to the maximum region. 

It is important to point out that the shower profile of the total energy release is 
not influenced when varying the simulation energy threshold within reasonable 
limits, as also shown in [22]. Adopting, for instance, higher thresholds would 
increase the contribution of the discarded particles, but this would just be 
balanced by a decreased contribution of the explicitely tracked particles. 



4 Lateral spread and particle densities 

For laboratory measurements of the fluorescence yield, the question arises 
whether the particle densities in air showers are small enough to allow for an 
undisturbed de-excitation or whether an ionized or excited air molecule might 
be influenced by another nearby shower electron. For this reason, the lateral 
spread of the energy release and the corresponding density of electrons and 
positrons as the component which dominates the energy release, is studied. 
More detailed investigations on the lateral distribution of the energy release in 
showers including the possibility of observing the shower width with fluores- 
cence telescopes are given elsewhere [30]; here we focus on questions relevant 
to laboratory measurements of the fluorescence yield. 
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Fig. 2. Upper panel: Energy release contained within a given distance of the shower 
axis versus the distance at shower ages s = 0.7 and s = 1.0. Lower panel: Lateral 
distribution of electrons plus positrons. 

In Figure 2 the amount of released energy contained within a given distance 
to the shower axis is plotted as a function of the distance (upper panel) for 
different values of the shower age s, defined as s = 3X/{X + 2Xmax)- Also 
given is the density in particles per square-meter of electrons and positrons 
(lower panel). As the lifetime of the radiating molecular states is of the order 
30—70 ns and thus comparable to the traversal time of the shower for a given 
air particle (or correspondingly to the shower "thickness" ) , transverse particle 
distances are regarded for a conservative estimate of the particle separation. 

Particles with core distances below 1 m contribute only little to the energy 
release: Though the particle densities are largest here, due to phase space 
the absolute particle number is comparatively small. The main energy release 
of about 80—85% occurs at distances between 1 — 100 m from the shower 
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axis. For core distances around 20 m, a typical transverse particle separation 
of ~ 1 mm for 10 EeV showers is obtained. This value holds for the shower 
maximum, at other stages of the shower development observable in the fioures- 
cence hght the particle separations are larger. The densities are roughly scahng 
with the primary energy, while the dependence on the primary particle type 
is of minor importance in this context. For 100 EeV showers at relatively 
small core distances of about 3 m, for instance, the average transverse par- 
ticle separation amounts to ~ 0.1 mm. Thus, with respect to the size of the 
ionization region around the charged particles in air (usually <1 /xm), this is 
a large separation resulting in a relatively "undisturbed" de-excitation of the 
air molecules. High-density particle bunches should therefore be avoided in 
fluorescence yield measurements as the fluorescence yield, apart from extreme 
cases of high shower particle densities, would be obtained in conditions not 
typical for fluorescence light emission in air showers. 



5 Energy spectra 

The particle energies contributing to the energy release, an important input 
quantity for the layout of fluorescence yield measurements, arc shown in Fig- 
ure 3. For this graph, only electrons and positrons are taken into account, 
since they dominate the energy release. The main contribution comes from 
particles with energies below 1 GeV, with a tail towards small energies. A 
broad maximum is visible at particle energies around 10—50 MeV, below the 
critical energy of electrons in air (~84 MeV). 

The spectral shape mainly reflects the particle energy spectrum [21]. Espe- 
cially the contributions of the lower energies are more pronounced, however. 
This is due firstly to the increased specific energy loss (Bethe-Bloch formula), 
and secondly to a larger average path length dx = dX/ cos{d) through the 
slant depth interval dX, since at lower electron energies the dispersion of par- 
ticle angles is increasing. As a guideline, in Table 1 the contribution to the 
electromagnetic energy release around shower maximum is estimated for dif- 
ferent ranges in kinetic energy. The value for E^in < 0.1 MeV is given by the 
contribution of the particles discarded due to the simulation energy threshold. 

Table 1 

Estimates for the contribution of different ranges in kinetic energy to the energy 
release by electrons and positrons. The uncertainty of the values is about ±2 (in %). 



Energy in MeV 


< 0.1 


0.1-1 


1-10 


10-100 


100-1000 


>1000 


Contribution in % 


10 


12 


23 


35 


17 


3 
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Fig. 3. Contribution to the energy release per matter traversed in shower direction 
as a function of the kinetic particle energy. Simulation for primary iron, 10^^ eV. 
Upper panel: Individual and combined contributions of electrons and positrons at 
shower maximum. Lower panel: Combined contribution of electrons and positrons 
for different shower ages (normalized to the same height of maximum). 



While at higher kinetic energies (-Efein > 300 MeV) electrons and positrons 
contribute about equally to the energy release, at lower energies only electrons 
survive due to positron annihilation. The annihilation photons will eventually 
transfer the energy to electrons by Compton scattering. 

The range of particle energies that mainly contribute to the energy release, is 
to a good approximation quite independent of the primary particle type (in- 
cluding primary photons) and primary energy. Also the dependence on shower 
age is small, as shown in the lower panel of Figure 3. For instance, at earlier de- 
velopment stages before the shower maximum, the contribution is only slightly 
shifted to higher electron energies. These results may be understood, since the 
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particle energy spectrum is known to show a small, but in this context only 
minor dependence on the primary particle type and on shower age [21,28]. 



6 Conclusion 

The presented method of treating sub-threshold particles allows the precise 
calculation of the total energy deposit. For shower calculations, the energy 
release provided by CORSIKA can be transformed to fluorescence light based 
on existing and upcoming fluorescence yield measurements. 

The energy release in air showers has also been studied with respect to cur- 
rently planned laboratory fluorescence yield measurements. Most relevant is 
the determination of the yield for electrons and positrons with energies in the 
range from sub-MeV up to a few hundred MeV. The typical lateral particle 
separation is relatively large with 1 mm or more for 10 EeV showers at shower 
distances where most energy is released and thus, presumably, the main fluo- 
rescence light production occurs. 

Finally, we want to mention that based on the concept of total shower size 

one can define the mean energy deposit per particle [1,7]. This quantity has 
been subject of many discussions and will be studied in detail in a forthcoming 
publication. 
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